Emission times and opacities from interferometry 
in non-central Relativistic Nuclear Collisions 
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The nuclear overlap zone in non-central relativistic heavy ion collisions is azimuthally very 
asymmetric. By varying the angle between the axes of deformation and the transverse direction of 
the pair momenta, the transverse HBT radii oscillate in a characteristic way. It is shown that these 
oscillations allow determination of source sizes, deformations as well as the opacity and duration of 
emission of the source created in any non-central high energy nuclear collisions. The behavior of 
the physical quantities with centrality of the collisions is discussed — in particular changes caused 
by a possible phase transition to a quark-gluon plasma. 

PACS numbers: 25.75 



Understanding particle emission in relativistic nuclear 
collision is crucial in order to determine whether a phase 
transition to quark-gluon plasma has occurred. Particle 
interferometry was invented by Hanbury-Brown & Twiss 
(HBT) for stellar size determination H and is now em- 
ployed in nuclear collisions J2-@|. It is a very powerful 
method to determine the 3-dimensional source sizes, life- 
times, duration of emission, flow, etc. of pions, kaons, 
etc. at freeze-out. Since the number of pairs grow with 
the multiplicity per event squared the HBT method will 
become even better at RHIC and LHC colliders where 
the multiplicity will be even higher. In this paper an ex- 
tension of the HBT method is presented for asymmetric 
and opaque sources created in non-central collisions. It 
is shown that combining HBT with a determination of 
the reaction plane can be exploited to find not only the 
size and deformation of the source but also its opacity as 
well as the duration of emission separately. 

The reaction plane, which breaks azimuthal symmetry, 
has been successfully determined in non-central heavy ion 
collisions from intermediate Q up to relativistic energies 
pp|pl| . The particle spectra are expanded in harmonics 
of the azimuthal angle 4> event-by-event ]9|-p^] 
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where <pn is the azimuthal angle of the reaction plane. 
Assuming that experimental uncertainties in event plane 
reconstruction can be corrected for, each event can be ro- 
tated such that <f>R ■ — 0. The asymmetry decreases with 
centrality (see Fig. 1) and vanish for the very central 
collisions, which are cylindrical symmetric. The stan- 
dard geometry has the z-direction along the longitudinal 
or beam axis and the x-direction along the impact pa- 
rameter b. Thus (x, z) constitutes the reaction plane and 
(x, y) the transverse direction with y perpendicular to the 
reaction plane. In HBT analyses one often parametrize 



the transverse particle distribution as a gaussian func- 
tion of transverse radius. It is therefore convenient for 
asymmetric sources to employ gaussian parametrizations 
in both transverse directions 



Sx(x,y) ~ exp 
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Here R x and R y are the gaussian transverse sizes of the 
collision zones at freeze-out. This elliptic source neglects 
directed flow and we shall therefore be restricted to re- 
gions around cms midrapidity where v\ vanish |lO|JTl| ]. 

The azimuthal asymmetry or "deformation" of the 
source can be defined as the relative difference between 
the gaussian radii squared 
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A simple estimate can be obtained from the full trans- 
verse extent of the initial overlap of two nuclei with radius 
Ra colliding with impact parameter b (see Fig. 1). They 
are simply: R x = Ra — b/2 and R y = \J R\ — b 2 /4, and 
the corresponding deformation is 
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The rms radii of the nuclear overlap zone weighted with 
longitudinal thicknesses results in a deformation that 
is only slightly smaller at semicentral collisions. How- 
ever, as the source expands the deformation S ~ (R y — 
R x )/(R y + R x ) decreases for two reasons. Firstly, the 
expansion increase {R x + Ry) and secondly, (R y — R x ) 
decrease because the average velocities are larger in the 
x- than y-direction. The latter is a consequence of the 
experimentally measured positive elliptic flow (v2 > in 
Eq. (|l|)) in relativistic nuclear collisions where shadowing 
is minor Jl(],[ll| . Measuring the decrease of the deforma- 
tion with centrality will reveal important information on 
the expansion up to freeze-out. For very peripheral col- 
lisions, where only a single nucleon-nucleon collision oc- 
curs, the source must be azimuthally symmetric, i.e., the 
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deformation must vanish and therefore Eq. (4), which 
assumes continuous densities, breaks down. 

The standard HBT method for calculating the Bose- 
Einstein correlation function from the interference of two 
identical particles is now briefly discussed. For a source 
of size R wc consider two particles emitted a distance 
~ R apart with relative momentum q = (pi — P2) and 
average momentum, P = (pi +p2)/2. Typical heavy ion 
sources in nuclear collisions are of size R ~ 5 fin, so that 
interference occurs predominantly when q^fi/R ~ 40 
MeV/c. Since typical particle momenta are p^P ~ 300 
MeV/c, the interfering particles travel almost parallel, 
i.e., pi ~ P2 — P 3> q. The correlation function due to 
Bose-Einstein interference of identical spin zero bosons 
7r ± 7r ± , K^K^, etc.) from an incoherent source is (see, 
e.g. 
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where S(x, P) is the source distribution function describ- 
ing the phase space density of the emitting source. 

Experimentally the correlation functions are often 
parametrized by the gaussian form 
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Here, q = pi — P2 = (q s , 9o, qi) is the relative momen- 
tum between the two particles and Ri,i = s,o,l,os,ol, si 
the corresponding sideward, outward, longitudinal, out- 
sideward, out-long and sideward-longitudinal HBT radii 
respectively. We have suppressed the P dependence. We 
will employ the standard geometry, where the longitudi- 
nal direction is along the beam axis, the outward direction 
is along Pj_ ~ p_i_,i, and the sideward axis is perpendic- 
ular to these. Usually, each pair of particles is lorentz 
boosted longitudinal to the system where their rapidity 
vanishes, y = 0. Their average momentum P is then per- 
pendicular to the beam axis and is chosen as the outward 
direction. In this system the pair velocity (3p—P/Ep 
points in the outward direction with (3 = pj_/mj_, where 
m± = \/ m 2 + p 2 ^ is the transverse mass, and both the 
R i and R s i vanish at midrapidity (see Ref. || for further 
analyses). Also R os vanishes for a cylindrically symmet- 
ric source or if the azimuthal angle of the reaction plane 
is not determined and therefore averaged over — as has 
been the case experimentally so far. The reduction factor 
A in Eq. (||) may be due to long lived resonances (||3), 
coherence effects, incorrect Coulomb corrections or other 
effects. It is A ~ 0.5 for pions and A ~ 0.9 for kaons. 

The Bose-Einstein correlation function can now be cal- 
culated for a deformed source. Let us first investigate 
transparent sources and, as in , parametrize the trans- 
verse and temporal extent by gaussians 



S{x,P) ~ S ± (x,y)exp[- 
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with longitudinal Bjorken flow, u — (cosh?y, 0, 0, sinhry). 
Effects of transverse flow will be discussed below. The 
transverse radii R x , R y are the gaussian radii at freeze- 
out, Tf is the freeze-out time and St the duration of emis- 
sion. 

In order to calculate the correlation function of (H) 
the gaussian approximation is employed (see, e.g., (J]) 
which results in a correlation function on the form as in 
Eq. (||). Inserting the source (0) in Eq. (||) and Fourier 
transforming we obtain the correlation function. Com- 
paring to the experimental parametrization of Eq. (^) , one 
calculates the HBT radii. For transparent sources the az- 
imuthal dependence of the HBT radii has been calculated 
in detail by Wiedemann JL4| . In the longitudinal center- 
of-mass system of the pair (y — 0), the HBT radii are 
around midrapidities 
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where R 2 — (R 2 + R 2 )/2 is the average of the source 
radii squared. As in asymmetric flow, Eq. ([!]), is the 
azimuthal angle between the transverse momentum p± 
and the reaction plane. It is therefore the angle by which 
the i? , s axes are rotated with respect to the (x, y) re- 
action plane (see Fig. 1). Near target and projectile 
rapidities the directed flow is appreciable and leads to 
cos0 terms in Eqs. (||||) The out- and sideward 

HBT radii show a characteristic modulation as function 
of azimuthal angle with amplitude of same magnitude 
but opposite sign. Measuring the amplitude modulation 
of R s ,o,os determines five quantities and thus overdeter- 
mines the three source parameters which are the source 
size R, deformation S and and duration of emission St. 

Next we consider opaque sources. In relativistic heavy 
ion collisions source sizes and densities are large and one 
would expect rescatterings. As a result particles are pre- 
dominantly emitted near the surface and arrive from the 
(front) side of the source facing towards the detector. In 
||| it was found that for opaque sources, where mean free 
paths are smaller than source sizes, X m fp^R, the side- 
ward HBT radius increase whereas the outward is signifi- 
cantly reduced. The simple geometrical cause is that the 
emission region in the outward direction is the surface re- 
gion which considerably narrower than the whole source. 
As in Glauber absorption is introduced by adding 

an absorption factor exp(— J o~p(x')dx') where a is the 
interaction cross section, p the density of scatterers and 
the integral runs along the particle trajectory from source 
point x to the detector. Defining the mean free path as 
A m /p = (cr/o(0)) _1 , where p(0) is the central density, the 
source is opaque when A m / P <C R and transparent when 
A m /p S> R. Calculating the correlation function for an 
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opaque source from Eq. @ and comparing to the defi- 
nition of the HBT radii in Eq. (||), one generally obtains 
by expanding for small deformations, 



R 2 a =g s R 2 [1 + 5 co$(2<f>)] , 
R? = 9o R 2 [1 - (5cos(20)] + (3 2 St 2 
R 2 os =g os R 2 5sm{2$). 
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Here g ,s,os are model dependent factors that are func- 
tions of opacity but independent of the deformation. For 
a gaussian source (p oc S±), which is moderately opaque 
(^mfp/R = 1), a numerical calculation gives g s ~ 1.4 
and g a ~ 0.9 (see also |l|]). For a disk source with trans- 
verse radius twice the gaussian radius 2R, that emits like 
a black body (i.e., X m fp *C R), one finds g s — 4/3 and 



4 (4-(f) ) -°-2 |15|. In all cases g os ~ g s . Gener- 



ally, (gs — g ) > and the difference increases with opac- 
ity. Only for a completely transparent sources is g s = g Q . 
In Fig. 2 the HBT radii of Eqs. (§|lj) are shown for a 
near-central collisions {5 — 0.2) with a moderate dura- 
tion of emission (f3 2 Sr 2 / g s R 2 = 1/4) for various opacities 
Knf P /R = 0.1, 0.5, 1.0, 2.0, oo. As the opacity increases, 
g /gs decreases and therefore also the outward HBT ra- 
dius and its amplitude. 

Comparing the HBT radii from an opaque source 
Eqs.( [l2f]l4| ) with those of a transparent source Eqs.(||- 
|l0| ), one notices that the amplitudes in R s and R differ 
by the amount (g s — g D ). The modulation of the HBT 
radii with (j) provides five measurable quantities which 
over-determinates the four physical quantities describing 
the source: its size R, deformation 5, opacity (g a — g s ) 
and duration of emission St, at each impact parameter. 
The azimuthal dependence of the HBT radii thus offers 
an unique way to determine the opacity of the source as 
well as the duration of emission separately. 

Experimentally, HBT analyses have not been combined 
with determination of the reaction plane yet. Conse- 
quently, the azimuthal angle <f> is averaged and informa- 
tion on three of five measurable quantities in Eqs. (|i2fjl4| ) 
is lost. From the angular averaged difference between the 
out- and sideward HBT radii 



(R 2 - R 2 )^ = /3 2 St 2 - (g s - 9o )R 2 
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one can only determine the sum of the positive duration 
of emission and negative opacity effect. Experimentally 
the difference is small; NA49 @ and NA44 [|| data even 
differ on the sign. Detailed analyses of the p± depen- 
dence of the HBT radii from NA49 data within opaque 
sources 16 indicate that the sources are transparent or 
at most moderately opaque. However, the NA44 data for 
which R ^R S necessarily requires an opaque source as 
seen from Eq. (|l5|). Furthermore, the p± dependence of 
the transverse HBT radii change if the source sizes, opac- 
ities, and duration of emission also are p_L-dependent. 

Transverse flow may affect the out- and sideward 
HBT radii as opacity, i.e., the factors g 0yS may de- 
pend on both. For transparent sources transverse flow 



has been studied in 14|i|Jl6| through the flow u = 
(7 cosh(?7), 7 sinh(^), uj_), where 7 = y/ 1 + u 'j_ , assuming 
that the transverse flow scales with transverse distance, 
~ uor±/R. Otherwise the same transparent gaussian 
source as in Eq. (0) was employed. To lowest order in 
the transverse flow both transverse HBT radii decrease 
by the same factor ~ (1 + u^m^/T) to leading order in 
Uq. This transverse flow correction is independent of the 
source size and therefore also the deformation. Conse- 
quently, spatial deformations reduce the amplitudes by 
the same amount in this model for both transparent (see 
fuj ) and opaque sources. There are, however, box shape 
models where the transverse flow reduce R more than 
R s fTq| . A similar flow effect is found in hydrodynamic 
models also jlTj although it is considerably less than the 
opacity effect, i.e., when the Cooper-Frye freeze-out con- 
dition is modified by removing the back side of the source. 
Preferrably, the p± dependence of the HBT radii should 
be measured as well as the transverse flow from appar- 
ent temperatures in order to determine the magnitude 
of the flow and opacity separately. The transverse flow 
might also be azimuthally dependent, i.e., gs,o 

depend 

on <f> and thereby change the amplitude. This effect can 
be estimated from the measured elliptic flow and we find 
that for semicentral Pb + Pb collisions at SPS energies 
( u x) ~ ( u y) ~ v 2> which is only a few percent. Therefore, 
the azimuthally dependence of the flow and its effect on 
the amplitude in the HBT radii is also of that order only 
which is much less than 5 ~ 0.5 for semi-central colli- 
sions. The conclusion is that besides opacity also trans- 
verse flow may affect the factors <? 0jS - but independently 
of azimuthal angle. Therefore Eqs. ( |l2f|i4| ) are still valid 
and can be used to extract the duration of emission un- 
ambiguously. 

Measuring the centrality or impact parameter depen- 
dence of the source sizes, deformation, opacity, emission 
times and duration of emission is very important for de- 
termining how the source change with initial energy den- 
sity. If no phase transition takes place one would expect 
that source sizes and emission times increase gradually 
with centrality whereas the deformation decrease approx- 
imately as (|3j). In peripheral collisions source sizes and 
densities are small and few rescatterings occur. There- 
fore, the source is transparent and the HBT radii are 
given by Eqs. (§|ll]). For near central collisions sources 
sizes and densities are higher which leads to more rescat- 
terings. Thus the source is more opaque and the ampli- 
tudes should differ. It would be interesting to observe 
this gradual change in amplitude with centrality. At the 
same time it would provide a direct experimental deter- 
mination whether the source is transparent or opaque 
as well as extracting the magnitudes of the opacity and 
duration of emission. 

If a phase transition occurs at some centrality, where 
energy densities exceed the critical value, one may ob- 
serve sudden changes in these quantities. The emission 
time and duration of emission has in hydrodynamic cal- 
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dilations jl7j been predicted to increase drastically lead- [4 
ing to very large Ri and R . A long lived mixed phase 
would also emit particles as a black body and thus the [5 
opacity should also become large. If droplet formation 
occurs leading to rapidity fluctuations, one may be able [6 
to trigger on such fluctuations and find smaller longi- 
tudinal and sideward HBT radii Furthermore, if P. 
an interesting change in these quantities should occur at 
some centrality, it would also be most interesting to look 
for simultaneous J /\& suppression, strangeness enhance- 
ment, decrease in directed, elliptic or transverse flow p9| , 
or other signals from forming a quark-gluon plasma. 

In summary, the importance of measuring the reaction 
plane and HBT radii simultaneously has been stressed. 
The modulation of the HBT radii with azimuthal angle 
between the reaction plane and particle transverse mo- 
menta can be exploited to obtain source sizes, deforma- 
tions, life-times, duration of emission and opacities sepa- 
rately. Tracking these physical quantities with centrality 
will provide detailed information about the source cre- 
ated in relativistic nuclear collisions and may reveal the [14 
phase transition to quark-gluon plasma. [15 
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FIG. 1. Reaction plane of semi-central Pb + Pb collision for impact parameter b = Rpb — 7fm. The overlap zone is deformed 
with R x < R y . The reaction plane (x, z) is rotated by the angle <j> with respect to the transverse particle momentum p± which 
defines the outward direction in HBT analyses. 
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FIG. 2. HBT radii vs. angle between reaction plane and transverse particle momenta. The HBT radii are normalized 
to the angle averaged sideward HBT radius squared, g s R 2 , for slightly deformed S — 0.2 source with duration of emission 
6%8t 2 I g s R 2 = 1/4 and with various opacities. The sideward HBT radius (full curve) is then the same for both transparent (Eq. 
(g)) and opaque (Eq. (^)) sources and likewise for the out-side HBT radii (Eqs. (Rl and (|l4[), dashed curve). The outward 
HBT radii are shown with chain-dashed curves for a gaussian source (Eq. (h3h and(foh for various opacities (from below and 
up): X m fp/R = 0.1, 0.5, 1.0, 2.0, oo. 
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